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Abstract 

A class of shape-invariant bound-state problems which represent transition 
in a two-level system introduced earlier are generalized to include arbitrary 
energy splittings between the two levels as well as intensity-dependent interac- 
tions. We show that the couple-channel Hamiltonians obtained correspond to 
the generalizations of the nonresonant and intensity-dependent nonresonant 
Jaynes-Cummings Hamiltonians, widely used in quantized theories of laser. 
In this general context, we determine the eigenstates, eigenvalues, the time 
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evolution matrix and the population inversion matrix factor. 
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I. INTRODUCTION 



The integrability condition called shape-invariance originates in supersymmetric quan- 
tum mechanics fl],@]. The separable positive-definite Hamiltonian Hi = A^A is called shape- 
invariant if the condition 



A{a 1 )A\a 1 ) = A\a 2 )A(a 2 ) + fl(ai) 



(1.1) 



is satisfied ||. In this equation % and a 2 represent parameters of the Hamiltonian. The 
parameter a 2 is a function of a x and the remainder R(ai) is independent of the dynamical 
variables such as position and momentum. Even though not all exactly-solvable problems 
are shape- invariant [[|, shape invariance, especially in its algebraic formulation has 
proven to be a powerful technique to study exactly-solvable systems. 

In a previous paper || we used shape-invariance to calculate the energy eigenvalues and 
eigenf unctions for the Hamiltonian 



H = A^A + - [i,i f ' 



1^3 



1) + Vm(a + A + 



where 



2 (^i ± ^2) 



;i-2) 



1.3) 



and with z = 1, 2, and 3, are the Pauli matrices. 

This is a generalization of the Jaynes-Cummings Hamiltonian |§. A different, but related 
problem was considered in Ref . |1(J . Our goal in this paper is to study further generalizations 
of the Jaynes-Cummings Hamiltonian, first by introducing a term proportional to 03 with 
an arbitrary coefficient (the so-called nonresonant limit) and then by taking into account the 
dependence of the coupling on the intensity of the field (the so-called intensity-dependent 
nonresonant limit). In addition to the energy levels we study the time evolution and the 
population inversion matrix factor. 

Introducing the similarity transformation that replaces ai with a 2 in a given operator 



and the operators 



T(ai) 6(ai) f \a x ) = 6{a 2 ) 



B + = it(ai)T(ai) 



:i-4) 



;i.5) 



B_=B\=f\a l )A{a l ), 
the condition of Eq. (|1.1|) can be written as a commutator 

[B_,B + ] = f t (a 1 )ii;(ai)f(ai) = R(a ) 
where we used the identity 

R(a n ) = f(a 1 ) i?(a n _i) T f (ai) , 



;i.6) 



;i.7) 



;i.s) 
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valid for any n. The ground state of the Hamiltonian H\ = A^A 
condition 

A | ^o) = o = zL | Vo) ; 

and the unnormalized n-th excited state is given by 
with the eigenvalue 

£ n = ^ i2(o fc ) . 
fe=i 



B+B_ satisfies the 



;i.9) 



We note that the Hamiltonian of Eq. (|1.2| ) can also be written as 



H 



f 


" 




■ft 


" 





±1. 




. 


±1 



;i.io) 



;i.n) 



;i.i2) 



where 



h ± = B + B_ + -R(a ) (ct 3 + 1) ± Vhtl [&+B- + a_B + 



;i.i3) 



II. THE GENERALIZED NONRESONANT JAYNES-CUMMINGS 

HAMILTONIAN 

The standard Jaynes-Cummings model, normally used in quantum optics, idealizes the 
interaction of matter with electromagnetic radiation by a simple Hamiltonian of a two- 



level atom coupled to a single bosonic mode fll"|-4t~6|j . This Hamiltonian has a fundamental 
importance to the field of quantum optics and it is a central ingredient in the quantized 
description of any optical system involving the interaction between light and atoms. The 
Jaynes-Cummings Hamiltonian defines a molecule, a composite system formed from the cou- 
pling of a two-state system and a quantized harmonic oscillator. In this case, its nonresonant 
expression can be written as 



H 



i f i + ]- [A, i f ] (a 3 + 1) + a (a + A + <r_ A f ) +hAa 3 , 



(2.1) 



where a is a constant related with the coupling strength and A is a constant related with 
the detuning of the system. 

However, the harmonic oscillator systems, used in this context, is only the simplest 
example of supersymmetric and shape-invariant potential. Our goal here is to generalize 
that Hamiltonian for all supersymmetric and shape-invariant systems. With this purpose 
and following Ref. M we introduce the operator 



S = a + A 



(7. 



.it 



(2.2) 
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where the operators A and A* satisfy the shape invariance condition of Eq. (|1.1|) . Using this 
definition we can decompose the nonresonant Jaynes-Cummings Hamiltonian in the form 



H = H D + H 



where 



int i 



H = S 2 , 

H int = « S + HA a 3 . 



(2.3) 



(2.4a) 
(2.4b) 



First, we search for the eigenstates of S 2 . In this case it is more convenient to work with its 
.B-operator expression, which can be written as [H] 



(2.5) 



where H 2 = TB_B + T^ . Note the freedom of sign choice in Eq. ( |2.5| ), which results in two 
possible decompositions of S 2 . Next, we introduce the states 



"T 


" 




. 


±1. 





-B + " 




■ft 


" 




H 2 


" 


B + B_ _ 




. 


±1. 




. 


Hi. 



*C±)) 



T 
±1 



I n) 



(2.6) 



where = [R(a%), Rfa), #(03), • • •] are auxiliary coefficients and, | m) and | n) 

are the abbreviated notation for the states | ip m ) and | ip n ) of Eq. (|1.10| ). Using Eqs. Ql.7|), 
Q2.5|) and ( |2.6|) , the commutation between f/i and a function of R(a,k), and the T-operator 
unitary condition, we get 



S 2 I 



"T 


" 




.0 


±1_ 




"f 


" 




.0 


±1. 





B + B_ + R{a ) 1 TCW I m) 
J [ I ra) 

^ m + i?(ao) OirC^lm) 
^nJLcWjn) 



(2.7) 



And using Eqs. ( |1.8| ) and ( |1.11| ) we can write 

f [£ m + R(a )} f t = f [i?( ai ) + fl(a 2 ) + ■ • • + i?(a m ) + i?(a )] ft 

= R(a 2 ) + R(a 3 ) H h -R(a m+ i) + i?(ai) = £ m+1 



Hence the states 



m I 



T 
±1 



a 



Cg> I m) 

771 + 1) 



(±) 
m+1 



m = 0,1,2, 



are the normalized eigenstates of the operator S 2 



(2.8) 



(2.9) 



(2.10) 



We observe that the orthonormality of the wavefunctions imply in the following relations 
among the C's 



5 



a 



(±) 



771+1 



\ 717 I m / 777 1 

(^m ) I = w - ^iciri = . 



(2.11a) 
(2.11b) 



Since S 2 and Hj nt commute then it is possible to find a common set of eigenstates. We 
can use this fact to determine the eigenvalues of Hint and the relations among the C"s 
coefficients. For that we need to calculate 



H 



tut 



#(±)\ = A (±) I # (±) ) 

771 / v 777 I 771 / 



where are the eigenvalues to be determined. Using Eqs. ( [2~2] 
eigenvalue equation can be rewritten in a matrix form as 



(2.12) 

and (pT9|) , the last 



P TB_ 
5+ft -p 



T 
±1 



CW | m) 



tf£> I m) 
CiSi |m+l> 



(2.13) 



where /3 = HA/ a. Since the C's coefficients commute with the A or operators, then the 
last matrix equation permits to obtain the following equations 



a 
a 



P - A^] (fC^ft) f | m) ± a C^if 5_ | m + 1) = 

Cjg x |m + 1) = 0. 



Introducing the operator [0] 



±) 

m 



-1/2 - 



one can write the normalized eigenstate of iZi as 

| m) = (Q f ) m | 0) 
and, with Eqs. ( |2.15| ) and ( |2.16| ) we can show that [ 



B+\m) = JS m+ i | m + 1) , 



TB- | m + 1) = y£ m+ i T | m) . 
Substituting Eqs. G2.17Q into Eqs. (|2.14|) we have 



ap - AW] (TCWft) ± a^S^Ctli \f\m) = 



a 



£ m+1 (fcWp)T[aP + >£ ) 



CM™ + i) = o- 



From Eqs. ( 2.18 ) it follows that 



AW = ±aJ£ m+1 + /? 2 



(2.14a) 
(2.14b) 



(2.15) 

(2.16) 

(2.17a) 
(2.17b) 

(2.18a) 
(2.18b) 

(2.19) 



and 
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a 



(±) 

771+1 



£ m+1 + /3 2 + /3 



ro+1 



^2.20) 



Eqs. (|m| ) and ( ggUp imply that 



(2.21) 



and the eigenstates and eigenvalues of the generalized nonresonant Jaynes-Cummings Hamil- 
tonians can be written as 



= £ m+l ± a/« 2 £ m+ i + h 2 A 



[2.22) 



and 



*(±)\ 



T 
±1 



c£0 1 m ) 

CW I m+l) 



m = 0,1,2, 



(2.23) 



a) The Resonant Limit 

From these general results we can verify two important and simple limiting cases. The 
first one corresponds to the resonant situation, for which A = (/3 = 0). Using these 
conditions in Eqs. fl2~20l) and ( ^22|) and Eqs. (TOD we get 



^m+1 ; 



and 



(2.24) 



a 



(±) 

m+l 



f/7(±)Tt = rf(±) 

J- * — • n -> * — ' im 



v/2' 



Therefore the Jaynes-Cummings resonant eigenstate is given by 



m ' V2 



T 
±1 



| m) 
m+l) 



m = 0,1,2, 



(2.25) 



(2.26) 



These particular results are shown in the Ref . |§ . 
b) The Standard Jaynes-Cummings Limit 

The second important limit corresponds to the standard Jaynes-Cummings Hamiltonian, 
related with the harmonic oscillator system. In this limit we have that T = T' — ► 1, 
B_ — > a, B + — ► a', A = uj — uj and £ m+ i = (m + l)hu. Using these conditions in the 
Eqs. (ET2CD , O) and Eqs. (gUD we conclude that 



= (m + ± Ja 2 hu(m + 1) + ft 2 (u; - w ) 



(2.27) 



and 
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M±) = „(±)/7(±) = n{T) 

^m+l Ira ^ ra w m 



(2.28) 



where 



i in 



1 + 5^ =f <5 m , 



(2.29a) 
(2.29b) 



(m + l)a: 2 /iCo> 

Therefore the standard Jaynes-Cummings eigenstate, written in a matrix form, is given by 



m / 



1 



"1 " 




m) 


.0 ±7^ ±} - 




. | m + 1) . 



m = 0,l,2,... 



(2.30) 



These results are shown in many papers, in particular, in the Ref. |L7 



III. THE TIME EVOLUTION OF THE NONRESONANT SYSTEM 

To study the time-dependent Schrodinger equation for a Jaynes-Cummings system in 
nonresonant situation 







ih- | = (h + eu) I *(*)) 

we can write the wavefunction as 

| (t)) = exp (-ik t/fi) | , 



(3.1) 



(3.2) 



and, by substituting this into Schrodinger equation and taking into account the commutation 
property between H D and Hj nt , we obtain 



We introduce the evolution matrix Uj(t, 0): 

| = Ui(t,0) | *i(0)). 

which satisfies the equation 

ih^-Ui(t,0) = H int tii(t,0), 



that is, in matrix form, written as 

rr, & ^ 2 

^1 ^2 



a 



5+ft -13 







U 12 ' 






U22. 



(3.3) 



(3.4) 



(3.5) 



(3.6) 
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where the primes denote the time derivative. One fast way to diagonalize the evolution 
matrix differential equation is by differentiating Eq. ( |3.5|) with respect to time. We find 



q2 



ih 



HLu,(t,o) 



which can be written as 



TJH TJII 1 

u ll u 12 

fj" JJ" 

u 21 u 22 



CO! 

uo 2 



Uu U 12 
U21 U 22 



where 



Tubt = a\JTB_B + T^ + (3 2 = V a 2 H 2 + (HA) 



huj 2 = a\jB + B_ +f3 2 = yja 2 E x + (HA) 2 . 



(3.7) 



(3.8) 



(3.9a) 
(3.9b) 



Now, since by initial conditions Uj(0, 0) = I, then we can write the solution of the evolution 
matrix differential equation ( |3TT| ) as 



Ui(*,0) 



cos (u>it) sin (uit) C 
sin {Cj 2 t) D cos {oJ 2 t) 



and the C and D operators can be determined by the unitarity conditions 

ut(t,o)u i (t,o) = u i (t,o)ut(t,o) = i. 



In the appendix A we show that the unitarity conditions (|3.1I|) imply 



c 



(H 2 



1/4 



TB_ 



(3.10) 

(3.11) 

(3.12a) 
(3.12b) 



D = —C^ . 

Therefore, we can write the final expression of the time evolution matrix of the system as 



Ui(t,0) 



cos (edit) sin (u>it) C 
— sin {uj 2 t) cos (cb 2 t) 



(3.13) 



For Jaynes-Cummings systems an important physical quantity to see how the system 



under consideration evolves in time is the population inversion factor [[Q],|l3|,yl| , defined by 

W(t) = a+(t) a-(t) - a-(t) a+(t) = a 3 (t) , (3.14) 

where the time dependence of the operators is related with the Heisenberg picture. In this 
case, the time evolution of the population inversion factor will be given by 



dt 



1 

ih 



U](t,0) a 3 ,H Ui(t,0) 



(3.15) 



and since we have 
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^3,H 

then Eq. Q3.15 ) can be written as 



a 3 , S = —2a S cr 3 



— = — S(t)a 3 (t). 



(3.16) 



(3.17) 



We can obtain a differential equation with constant coefficients for cr 3 (t) by taking the time 
derivative of Eq. (|3.17 ) 



d 2 a 3 (t) 2ia \ dS(t) 



dt 2 



-T- -^^(t)+S(t) 



Having in mind that 



dS(t) 
dt 



ih 



dt 



Uj(t,0) 



(3.18) 



S,H 



Ui(t,0) 



and that 



we conclude that 



S,H 



a/? S,a 3 =2a/?S6"3 



dS(t) 2ia(3 



S{t)a 3 (t). 



dt h 

Using Eqs. (|3~T7| ) and (|3~2lD into Eq. (glgD we obtain 

d 2 a 3 {t) 



dt 2 



+ & 2 a 3 (t)=F(t) 



where 



& 2 = 4 Ss 2 



F(t) 



ft 

4a 2 (3 



Ut^OjSUi^O). 



(3.19) 
(3.20) 

(3.21) 

(3.22) 

(3.23a) 
(3.23b) 



Eq. ( |3.22D corresponds to a non- homogeneous linear differential equation for a 3 (t) with 
constant coefficients since S 2 and H commute and, therefore, is a constant of the motion. 
The general solution of this differential equation can be written as 

a 3 (t) = a H (t)+a p (t), (3.24) 

and each matrix element of the homogeneous solution, satisfies the differential equation 



dt 2 



+ P 2 af k (t) = 0, j, fc = l,or2 



(3.25) 



with 
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%v x = 2a\]TB_B + T^ = 2 V a 2 H 2 , (3.26a) 



hv 2 = 2a^B + B_ = 2V« 2 H x . (3.26b) 
The solution of Eq. ( |3.25| ) is given by 

^(*) = yi(*)cifc + %(*)4, ( 3 - 27 ) 

where 

Vj(t) = cos (C>jt) (3.28a) 

Zj(t) = sin {9 jt) , (3.28b) 

and the coefficients djk and djk can be determined by the initial conditions. 

The matrix elements of the particular solution of the 03 (t) differential equation need to 
satisfy 

d 2 a p (t) 

"r 1 + »j a fk(t) = F jk (t) , j, k — 1,ot2, (3.29) 



dt 

and they can be obtained by the variation of parameter or by Green function methods, 
giving 

= p j x [^m)hk{o - m) [ hkio) , (3.30) 

where we used that the Wronskian of the system of solutions Vj{t) and Zj(t) is given by Uj. 

After we determine the elements of the F(£)-matrix, it is necessary to resolve the integrals 
in Eq. ( |3.30|) to obtain the explicit expression of the particular solution. In the appendix 
B we show that, using Eqs. Q2.2Q , (|3.13|) , and ( |3.23|) , it is possible to conclude that these 
matrix elements can be written as 



° xx {t) = i^o x 1 JfB~{z 2 {t) G ( ckt^2^2^i)-Ht) 4steA^A)}#2 1/4 

+ i^P^H^ {z x (t)g^J(t;i> x ,u x ,u 2 ) -y x (t)g^(t;u x ,u x ,u} 2 )^ \] 5+Tt , (3.31a) 

= ^y/fS^fa® $cc(t; ^AA) - fa{t) Gsc{t; *> 2 AA)} \/fB~l 

+ ^r 1 ^ 74 {k(t) G { ss\t; z>! A, £2) + Vx{t) g { jj(t; P x , u x , Qj 2 )}Hl /A , (3.31b) 



= y> 2 l 4&4^ {M*) G { c + c(t; £i AA) - m{t) g ( sc(t; *>i AA)} 

+ i^2 X Hl /A [z 2 {t) Gss{t\ z> 2 , u 2 , ljx) - y 2 {t) Q^{t; z> 2 , tb 2 , lj x )} H 2 /A , (3.31c) 
aUt) = i^JIj^ {z x (t) Q™(t; z^AA) - y x {t) g { s + s \t; v x , u x , u 2 )} H{ /A 

+ i^u 2 l Hl /A [z 2 {t) g$(t; z> 2 AA) + Ht) gccfr i>2 A A)} y/fSl, (3.31d) 
where 7 = Aa 2 (3/fi 2 , and the auxiliary functions are given by 
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QxY^PA^r) = T X y{t\V- q,r) ± Fxyfa P + q, f) , X, Y = C or S , 



(3.32) 



with 



Fccit'iX, w) = / d£ cos(x^) cos(w)£) 
Jo 



^.2m^jj2n ^2m+2ri+l 



m % K } (2m)!(2n)!(2m + 2n + l) 

J~cs(t'i x,w) = / d£ cos(x£) sin(w£) 
Jo 

oo »2m ;2n+l j.2m+2n+2 

~ m jrJ 1 (2m)!(2n+l)!(2m + 2n + 2) 
J~sc{t) x,w) = / d£ sin (x£) cos (w£) 

oo -2m+l,;,2n i2m+2n+2 

" m j?J 1 (2m+l)!(2n)!(2m + 2n + 2) 



Fss{t; x,w)= f d£ sin (xg) sin («;£) 
Jo 



t 



2m+2n+3 



m>n=0 (2m+l)!(2n + l)!(2m + 2n + 3) 

With these results for the particular solution we can conclude that 

p, x dafJO) 



Now, using Eqs. (|3.17|) , (|3.24| ), ( |3.27|) , (|3.34j) and the initial conditions, we have 

[^(0)], 

'da 3 (0) 



lij - C ij 



dt 



J >j 



2ia 
IT 



S(0)<7 3 (0) 



V{ dij . 



(3.33a) 



(3.33b) 



(3.33c) 



(3.33d) 



(3.34) 

(3.35a) 
(3.35b) 



Therefore, the final expression for the elements of the population inversion matrix of the 
system can be written as 

[a-s^ly = cos (Pit) [a-3(0)] o . + ^ sin (u t t) zT 1 [s(0) <t 3 (0)] .. + af^t) . (3.36) 
Again, using these final results we can verify two important and simple limit cases, 
a) The Resonant Limit 



The first one corresponds to the resonant situation (A = 0). Eqs. (|3.9|) , (|3.13| ), (|3.26|) 
and ( p.31|) allow us to conclude that, in this case, the evolution matrix of the system is given 
by 



Ui(t,0) 



cos ( huit] sin ( i/>it ) C 



sin 



|^2M & COS (js^t 



(3.37) 
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and the elements of the population inversion of the system are 



[h{t)]ij = cos (Pit) [<x 3 (0)L, + ^ sin (utt) ur 1 [S(0) a 3 (0) 



(3.38) 



b) The Standard Jaynes-Cummings Limit 

This second important limit corresponds to the case of the harmonic oscillator system, 



and in this limit we have that T = T* 



1, B. 



a, B< 



V and [i 



a, a 



htu. 



With these conditions the operators (b\ and Cj 2 commute, and this fact permits to evaluate 
the integrals related with the particular solution of the population inversion elements using 
trigonometric product relations. Using that and the expressions obtained in the appendix 
B, after a considerable amount of algebra and trigonometric product relations we can show 
that is possible to write the expressions for the <r^(t)-matrix elements as 



^u(*) = 1 v / a{/C 5 (t; ^2,^1,^2) -K. 8 (t\ w 2 ,-Wi,/> 2 )} (aa f ) 



1/4 



- z^z>f 1 (aa f ) 1 {/Cs(t;u> 2 ,a>i,z>i) - JC s (t;u 2 , -a>i,i>i)} Vat 
= ^f 1 v / S{/C c (t; ^2,^1,^2) -£c(*;w 2 ,-&i,i>2)}Va 

- I^r 1 (aa f ) 7 {£<?(*; w 2 ,&i,i>i) -/C c (t; a> 2 ,-cDi,i>i)} (a f a) 7 



= 0^2 1 v / at{/C c -(t; ^ 2 ,wi,z>i) + /C c (t;a> 2 ,-^i,i > i)} Vat 



7.-1 /~t-V /4 



-p 2 iaa 



{K c {t; ui 2 , &i, h) - £c{t; u> 2l -u>i, h)} (atf 



1/4 



^22^) = «?^2 lv/ ^ {^s(*5 w 2 , wi, *>i) + /C s (*; lj 2 , -&\, h)} (a j aj 



1/4 



- z^z> 2 1 (a f a) 1/4 {JC s {t; u 2 ,lji,i> 2 ) + /C s (t; tD 2 , ^2)} v 7 ^ , 

where, now, the auxiliary functions are given by 

f sin [(p + q) t] — (p + g) sin (ft) 



JC s (t;p,q,r) 



K. c {t]p,q,f) 



r cos 



f 2 — (p + (?) 
■f g) t] — f cos (ft) 



f 2 — (p + qY 



(3.39a) 
(3.39b) 
(3.39c) 
(3.39d) 

(3.40a) 
(3.40b) 



Considering the expressions above we may easily verify that the particular solution for the 
population inversion factor must still satisfy the initial conditions fl3.34|) . Therefore, in this 
case the final expression for the population inversion factor has the same form given by 
Eq. ( ET36D , with 



Hui = aJaa^ + (3 2 



%u) 2 



aJ at a + f3 2 . 



(3.41a) 
(3.41b) 
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IV. THE GENERALIZED INTENSITY-DEPENDENT NONRESONANT 
JAYNES-CUMMINGS HAMILTONIAN 



A variant of the Jaynes-Cummings model takes the coupling between matter and the 



radiation to depend on the intensity of the electromagnetic field ||13|,|15|,|16yi8| . This model 
has great relevance since this kind of interaction means effectively that the coupling is pro- 
portional to the amplitude of the field which is a very simple case of a nonlinear interaction 
corresponding to a more realistic physical situation. The results of this model can also give 
insight into the behavior of other quantum systems with strong nonlinear interactions. In 
this section we generalize the standard intensity-dependent nonresonant Jaynes-Cummings 
model to a shape-invariant one. 

The expression for the intensity-dependent nonresonant Jaynes-Cummings Hamiltonian 
can be written as 



H = i f i+^ (<r 3 + 1) + a ( a+A \JA^A +a^A^A A^)+hA& 3 . 



(4.1) 



Note that here the constant a is dimensionless. To generalize Hamiltonian ( |4.1| ) for all 
super symmetric and shape-invariant systems, we can use the operator S, given by Eq. Q2.2| ), 
and further by introducing the following operator 



A \JAU +a^A^AA ] . 



(4.2) 



Again, the operators A and A^ satisfy the shape invariance condition, Eq. ( |1 . 1| ) . Using 
operators S and Sj we can decompose the Jaynes-Cummings Hamiltonian in the form 



H = H + H 



int i 



(4.3) 



where 



H 



int 



a Si + hA <5"3 . 



(4.4a) 
(4.4b) 



In this case, H int can be written as 



Hint = « 



P 



B+B. 



a 



\IHi 

-a 



(4.5) 



Here we can follow the same development of the section II, with the same notation. So, 
using Eqs. (pT2|) , (|2~9|), fl4~2] ) and (JO]), the eigenvalue equation 



can be written in a matrix form as 



a 



p TB_ \JH X 

'H x B+f f -P 



T 
±1 



a 



Cg> | m) 
mli I m + !) 



CW | m) 



(4.6) 



• (4-7) 
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Again, since the C's coefficients commute with the A or A 1 operators, then the last matrix 
equation permits to obtain the following equations 



a; 



(3 - A W (TC^T 1 ") f\m)±at C^TB. \j H x \ m + 1) = 



a (f C&T f ) \]H l B + \m)T [a/3 + 
Now, using from Eqs. Q2.15|) to (|2.17| ) we have 



a 



(±) 

m+l 



m + 1) = . 



'#i | m + l) = T5_ J£ m+1 | m+l) 



£ m+ i T5_ | m + 1) 

£ m +l ^ I m ) ^ 



and 



H 1 B + \ m) = y Hi J£ m +i | m + 1) 



£ m+1 yHi | m + l) 
£ m+ i | m + 1) , 



Using Eqs. (gj) and flEIfl|) , then Eqs. take the form 

{ [a/3 - AgO] (f C^f t) ± a S m+ i C^Jf | m) = 
{a £ m+1 (TCg* P) T [a/3 + C^} | m + 1) = 

From Eqs. (|4.11|) it follows that 



X^ = ± a j£l + i + P 



and 



a 



^+i + /5 2 + /3 



m+l 



m+l 



(TCg=)f t) . 



Eqs. and ( ^TT5|) imply that 



c {±) =( 7(+) 

^m+l w m 



(4.8a) 
(4.8b) 



(4.9) 



(4.10) 

(4.11a) 
(4.11b) 

(4.12) 
(4.13) 
(4.14) 



and the eigenstates and eigenvalues of the generalized intensity-dependent nonresonant 
Jaynes-Cummings Hamiltonian can be written as 



£g=> = £ m +l ± V« 2 £1+1 + 



and 



T 
±1 



C£0 | m) 
| m+l) 



m = 0,1,2, 



(4.15) 



(4.16) 
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a) The Intensity-Dependent Resonant Limit 



From these general results we can again verify the two simple limiting cases. The first 
one, corresponding to the resonant situation, is for A = (j3 = 0). Using these conditions 
into Eqs. ( |4.13|) and ( |4.15|) and Eqs. (|2.11|) we can promptly conclude that 



'l±a) £ 



m+l j 



(4.17) 



and 



a 



(±) 

m+l 



v/2' 



(4.18) 



Therefore the intensity-dependent resonant Jaynes-Cummings eigenstate is given by 



m ' s/2 



T 
±1 



| m) 
m + l) 



m 



0,1,2, 



(4.19) 



If we compare this last particular result with that one found in the reference |§, we con- 
clude that the intensity-dependent and intensity-independent generalized Jaynes-Cummings 
Hamiltonians have the same eigenstates in the resonant situation. 

b) The Standard Intensity-Dependent Jaynes-Cummings Limit 



The second limit, corresponding to the standard intensity-dependent Jaynes-Cummings 
case, is related with the harmonic oscillator system. In this limit we have that T = — > 1, 
B_ — > a, B + — ► a', A = uj — u and £ m +± — (m + l)tku. Using these conditions in 
Eqs. ( |4.13| ), ( |4.15| ) and Eqs. (|2.11|) we can promptly conclude that 



and 



= (m + l)hu ± hJa 2 uo 2 {m + l) 2 + (u - uo ) 2 



(4.20) 



M±) = ~(±W±) = /7(t) 

w m+l lm w m w m 



1+ 7: 



(4.21) 



where 



im 



1 + $m =F $m ■ 



UJ — UJ n 



aui{m + 1) 



(4.22a) 
(4.22b) 



Therefore the standard intensity-dependent Jaynes-Cummings eigenstate, written in a ma- 
trix form, is given by 



1 



1 




(±) 

in 



| m) 
m + l) 



m = 0, 1,2,--- . 



(4.23) 
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V. TIME EVOLUTION OF THE INTENSITY-DEPENDENT NONRESONANT 

SYSTEM 



To resolve the time-dependent Schrodinger equation for intensity-dependent nonresonant 
Jaynes-Cummings systems: 



H„ + H 



hit 



we can again write the state | ^(i)) as it is given by Eq. flO] 
to ( |3.5|) , we can write the matrix equation 



*(*)> (5-1) 
Then using from Eqs. (P-3|) 



ih 



TJ[ X TJ\ 



12 



TP TP 

u 21 u 22 



a 



(3 TB_ \JHx 

Hi S+ft -/3 



U21 U22 



To diagonalize this evolution matrix differential equation we can differentiate Eq. (|3 
respect to time. After that, if we again use the same Eq. ( |3.5|) , we find 



H m ^U l (t,0) = lHLU,(t,0) 



which can be written as 



TP' TP' 

v \\ u 12 

U" TP' 

u 21 u 22 



£>i 
Cj 2 



ih 



U21 U22 



where 



a 2 Hi + (hA) 2 



7ub x = a\l(TB_B + T^) 2 + (3 2 = 
hu 2 = a^(B+B_) 2 + P 2 = ^a 2 H 2 + (hA) 2 . 



(5.2) 
with 

(5.3) 

(5.4) 

(5.5a) 
(5.5b) 



Now, using the initial conditions Uj(0,0) = i, we can write the solution of the evolution 
matrix differential equation ( |5.3|) as 



Ui(t,0) 



cos (&it) sin (uit) C 
sin {uj 2 t) D cos {Cj 2 t) 



(5.6) 



where the C and D operators can be determined by Eq. ( |3.11| ). Following the same steps 
used in the appendix A, we can conclude that these operators must have the form given 
by Eqs. ( |3.12| ). So, in this case the final expression of the time evolution matrix \Ji(t, 0) is 
given by Eq. ( |3.13|) as well. 

To obtain the population inversion factor we can again follow the steps from Eq. ( J3.14 ) 
to ( |3.30|) , but replacing the operator S by the operator Sj. Besides that we have 



Tiv x = 2a TB^B + f ] = 2a H 2 
hu 2 = 2a B + B_ = 2aH 1} 



(5.7a) 
(5.7b) 



instead Eqs. (|3.26| ). Here, we can again use the development shown in the appendix B, 
just replacing S by Sj, to obtain the explicit form of the matrix elements for the particular 
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solution of the population inversion factor, given by Eq. ( |3.30| ). So for a shape-invariant 
intensity-dependent nonresonant Jaynes-Cummings system, these matrix elements are given 

by 



i— v 



i 1 yTB^ {z 2 (t) ^2,^2,^1) - fa{t) Gss\t; £> 2 , u 2 , &i)} H. 



3/4 



2 

+ I — V " 

2 
2 

7.-1 £3/4 

2 



(5.8a) 



= tyTWfB- fat) g { c ^(t; 02,£>2M) ~ hit) g { g$(t; v 2 ,Lb 2 M)} \jH 2 fB_ 

+ {^(t) z> 1; £ 1} £ 2 ) + ft(t) g { cs \t; h^u 002)} Hi' 4 , (5.8b) 

b + t \h 2 [hit) g$(t; *>iA,£ 2 ) - yiW *l A,£ 2 )} V^+Tt 

+ J^ 1 ^! 74 {«2(*) ^ssC*; #2> ^2, wi) - 0cs (*5 £2, w 2 , &i)} # 2 /4 , (5.8c) 



1/4 



+ O-v^BL 1 ^ fat) g$(t; z> 2 , w 2 , £1) + m(t) $cc{t\ h, & 2 , ^H 2 TB_ . (5.8d) 



Yet the auxiliary functions, {?xy(£;p, g, f), are given by Eqs. ( |3.32| ) and ( |3.33| ). From 
Eqs. (|3.34j ) and (|3.35|) , we have for the elements of the population inversion matrix: 



2ia r* n 

[<7 3 (t)]o- = cos fo^L + — sin (utf) z/r 1 S,(0) <x 3 (0) . . + . (5.9) 



a) The Intensity-Dependent Resonant Limit 

In this limit we set (A = 0), so the evolution matrix of the system is given by 

cos(|z>itj sin ( \vit\ C 



UlM) " • cos(lz> 2 t) 
and the elements of the population inversion factor can be written as 

[& 3 {t)]ij = cos (£>it) [ct 3 (0)] o . + — sin(^) z/r 1 [S i (0)<r 3 (0) 



(5.10) 



(5.11) 



b) The Standard Intensity-Dependent Jaynes-Cummings Limit 

For the case of a harmonic oscillator system (T = — ► 1, B- — ► a, B^ 
[a, a/] = huj), we have for the <T^(t)-matrix elements the following expressions 



^u(t) = 1 v / S{/C 5 (t;^ 2 ,^i,z> 2 ) -/C 5 (t;^ 2 ,-^i,z> 2 )} (aa f ) 



3/4 



^1 1 (aa f ) 3/ {/C 5 (t;^ 2 ,cji,z>i) - /C 5 (t; u> 2 , -a>i, z>i)} v 7 ^ 7 



a* and 



(5.12a) 



= 1 v / ^{/C c (t; 0)2,^1, z> 2 ) - K.c(t; u 2l -ui, u 2 )} V aata 

- («a f ) 7 {K-c(t; u 2 , &x, v\) - fCc(t; a>2, -£1, £1)} (a f a) 1 (5.12b) 



(7 



" P (<) = -^Vataat {/C c (£; cj 2 , a> 1; z>i) + /Cc(t; w 2 , i>\)} vfiT 
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2 

■aT 1 (ala) 1 4 {/C c (^^2,^i, Aj) - ^Co(*;w 2 ,-a)i,z>2)} (aa 1 ") 374 (5.12c) 



2 



a. 
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= i^z> 2 Va+aat {/C 5 (t; c2> 2 , Cj\, V\) + /C s (t; a> 2 , h)} (a^aj / 

— (o^aj ^ {fcs(t] &2, v-i) + fcs(t] — z> 2 )} V aa^a , (5.12d) 

where the auxiliary functions, JCs(t',p, q,f) and JCc(t',P, q,f), are given by Eqs. fl3.4U| ). The 
final expression for the population inversion factor has the same form given by Eq. ( |5.9|) 
with 

hui = 2a aaJ , hv 2 = 2a a^a , (5.13a) 

%Cj x = a^(m) 2 + (3 2 , huj 2 = a^J(a^a) 2 + (3 2 . (5.13b) 



VI. CONCLUSIONS 

In this article we extended our earlier work || on bound-state problems which represent 
two-level systems. The corresponding coupled- channel Hamiltonians generalize the nonreso- 
nant and intensity-dependent nonresonant Jaynes-Cummings Hamiltonians. In the case of a 
nonresonant system, if we take the starting Hamiltonian to be the simplest shape-invariant 
system, namely the harmonic oscillator, our results reduce to those of the standard nonres- 
onant Jaynes-Cummings approach, which has been extensively used to model a two-level 
atom-single field mode interaction whose detuning it is not null. In addition we have studied 
time evolution and population inversion factor of the both kind of generalized systems. 

These models are not only interesting on their own account. Being exactly solvable 
coupled-channels problems they may help to assess the validity and accuracy of various 



approximate approaches to the coupled- channel problems ]19| . 
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Appendix A 

Here we give the steps used to obtain the specific form of the operators C and D. Using 
Eq. ( |3.10| ) into the unitary condition equation ( p. 11] ) actually we can show that the C and 



D operators need to satisfy the following six conditions 



C& = C ] C = 1 (la) 

DD* = D^D = 1 (lb) 

sin (u 2 t) = — sin (&\t) C (lc) 

D cos (pit) = - cos (u 2 t) C 4 . (Id) 



At this point we can use the following property 



T5_ Cj 2 = yTS_ ^a 2 B + B_ + (3 2 /h 
a 2 TB_B + B_ + TB^f3 2 /h 



a 2 TB_B + T^TB_ + (3 2 TB_/h 



a 2 TB_B + T^ + (5 2 /h \JTB 



Coi yTB_ . (2) 



Then, with this result we have 



'TB_ Cj\ = \/TB_ Cj 2 Cj 2 = u>! \JTB_ Cj 2 = Cj{ yTB_ , (3) 
and finally, by induction, we conclude that 



<TB_ Qj n 2 = Q\ \JTB_ . (4) 
In the same way, 



5+Tt lo x = \JB + T^ \]a 2 TB_B + T^ + (3 2 /h 

a 2 B + f^TB_B + f^ + B + f^ 2 /h 
a 2 B + B_B + f^ + (3 2 B + fyn 



a 2 B + B_+(3 2 /fi \JB + T1 



&2 \]B + T\ . (5) 



Then, with this result we have 



'B + f t Qj\ = yjB+fi Lbx Cj x = u 2 y B + f^ u x = Qj\ yB + f^ , (6) 
and finally, again by induction, we get 



5+Tt <2£ = QjI y B + T^ . (7) 
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Using the properties given by Eqs. (|]) and (|7|) and the forms of C, D operators, defined by 
Eqs. ( |3.12j ), we can verify that 



'-i) 1 f— 1 



H^ H£ 



and 



&c = bfr = V^yi im = 4^4= V^t^ = 1 • ( 9 ) 

Also using the series expansion of the trigonometric functions, we can show that 



D^&t) = w JTB-U-l) (2n+1)! 



2n+l 



E(- 1 ) n k^ruV^ 



4 1/4 n^o ' (2n + l)! 



OC 



Wit) — I 



= -wn(uiit)6, (10) 

where we used the commutation between H 2 and Ci}\ (see Appendix B). In the same way 
we can prove that 



D cos(^) = 4^-^, E(-l) n ^ 

(o>it) 



2/i, 



^ 2 1/4 ^o V ' (2n)! 

\2n 



S(-i)V^- (2jl)! ^ 

E(-i)^v^- 7 



(2n)! v + 

- cos (w 2 t) (7 t . (11) 



Again, we used the commutation between H 2 and Cj\. 



Appendix B 



In this appendix we show the necessary steps to obtain the explicit expressions of the 
particular solution elements of the population inversion factor. To resolve the integrals in 
Eq. (|3.30|) , first we need to determine the elements of the F(i)-matrix. To do that we can 
use Eqs. flTg) , ( gjgg) , and ( 0H to write down 
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F n (t) = -7 {cos (u x t) TB_ sin {uj 2 t) & + C sin {u 2 t) B + T^ cos (u>i£)j 

= 27 |\/t_B_ cos (uj 2 t) sin (c2>it) iY 2 ^ 4 — H 2 4 sin (cDxt) cos (<£ 2 i) \J B+T^ 
F\ 2 {t) = 7 {cos {(bit) TB_ cos (t2> 2 t) — C sin (u) 2 t) B + T^ sin (coit) C\ 

= 7 |\/t_B_ cos (u) 2 t) cos (o)i£) \ItB- + -£T 2 ^ 4 sin sin (u> 2 £) -ff 1//4 | 
F 21 (t) = 7 {cos (<j 2 t) 5 + T f cos pit) - (7 f sin (&!*) TB- sin (£ 2 £) C 1 } 

= 7 j\/-B + Tt cos (o>it) cos {u) 2 t) \/ B + T^ + i7^ 4 sin (a> 2 t) sin (c^i) i^ 1 / 4 
p 22 {t) = 7 {(7+ sin (uit) TB_ cos (w 2 t) + cos (u 2 t) B + f ] sin (wrf) d} 

— ij \ \J B + T^ cos (u>it) sin (o> 2 t) H\^ a — H^ 4 sin (a) 2 t) cos (u>ii) \JtB_ 



(12a) 



(12b) 



(12c) 



(12d) 



where 7 = Aa 2 (3/fi 2 . Here we used the properties ([I]), (0j) and (|7|), together with the 
following operators relations 



B+f t C=-& \JfB_ = iH{ /A 



(13a) 
(13b) 

Now, keeping in mind that %] = 0, (j = 1, or 2), so we may use the trigonometric rela- 
tionships involving products of trigonometric functions with arguments Ujt and Ujt (since we 
have exp (i>jt) exp (±%i) = exp \{0j ± Uj)t\). Then, using those relationships, the following 
commutators 







(14) 



and the same properties ([[]), (|J) and (^), we can show that 



■ 7 



2/i 00 = i^\jTB_ {cos [(z> 2 - cb 2 )t] sin (a>ii) + cos [(r> 2 + d> 2 )£] sin (o>it)} i7 2 



1/4 



+ i^-i^ 4 {sin [{pi — Co\)t] cos (cu 2 t) — sin \{p\ + c!>i)£] cos (u> 2 t)} y B + Ti (15a) 
Vi(t) F 12 (t) = ^\JfB^ {cos [{02 - Cj 2 )t] cos (a>it) + cos [{0 2 + u) 2 )t] cos (&it)} \JfB_ 

+ ^--ff 2 ^ 4 {sin [(z>! + u>i)£] sin {u 2 t) — sin [(j/j — &i)t] sin (a) 2 t)} i/ 1 ^ 4 (15b) 
fa{t) F 21 {t) = -y/B+fi {cos [{0! - Cj x )t] cos (w 2 t) + cos [{v x + &i)t] cos (w 2 t)} \/-B + ft 

+ ^i7 1/4 {sin [(z> 2 + & 2 )i] sin (a)it) - sin [(z> 2 - uj 2 )t] sin (o>it)} # 2 /4 (15c) 



• 7 



y 2 (t) F 22 (t) = i-y B + T^ {cos [(z>i - u>i)t] sin {u 2 t) + cos [(i>i + cD x )t] sin (u) 2 t)} 



1/4 



+ C^H.\^ {sin [(z/ 2 — u) 2 )t] cos (u>ii) — sin [{0 2 + a) 2 )t] cos {&it)} \JtB_ . (15d) 



In a similar way, we can show that 
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zi(t) F n (t) = i?U]TB- {sin [(z> 2 - ^a)*] sin (u^i) + sin [(z> 2 + £> 2 )£] sin ((Jit)} # 2 /4 



— z-^-i^ 4 {cos [(z>i — u>i)t] cos (u> 2 £) — cos [(z>i + Ct>i)t] cos (ca>2*) } V B + T^ (16a) 
21 (*) Aa(*) = ?V^^ { sin [(^2 - w 2 )t] cos pit) + sin [(z> 2 + w 2 )t] cos (u)it)}yJ¥5L 

— t^H^ 4 { cos [(Al + sin (u) 2 t) — cos [(Pi — (2>i)i] sin (a) 2 t) } (16b) 
z 2 (*) = ~^B+fi {sin [(r>i - Cj x )t] cos (w 2 t) + sin [(i>i + cDi)t] cos (o> 2 *)} ^B+ft 

— ^-^i^ 4 { cos [(^2 + &2}t\ sin (u)it) — cos [(z> 2 — u> 2 )t] sin (uit) } H 1 ^ (16c) 



2 2 (t) F 22 (t) = i— y 5 + Tt {sin [(i>\ - Ui)t] sin (u 2 t) + sin [(z>i + ui)t] sin (u) 2 t)} -f/y 



— i^H^ 4 {cos [(z> 2 — u) 2 )t] cos (d>it) — cos [(z> 2 + o) 2 )i] cos (£>it)} \JtB_ . (16d) 

The non-commutativity between the operators a)i and cD 2 imply that to calculate the 
integrals involving the terms given by Eqs. fljjD and (|16|) we need to use the series expansion 
of the trigonometric functions. In this case the integrals can be easily done because the time 
variable can be considered as a parameter factor. Finally, using these results into Eq. (|3.30| ) 
is trivial to find the expression (|3.31|) for the matrix elements of the particular solution. 
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